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SUMMARY 


The  theoretical  analysis  of  synthetic  aperture  techniques  presented  in  Memoran¬ 
dum  4532  is  applied  to  the  inverse  synthetic  aperture  radar  (ISAR)  imaging  of  ship 
targets.  The  point  spread  function  of  an  ISAR  system  is  derived  and  discussed  in 
detail.  The  explicit  representation  of  phase  information  within  our  formalism  makes 
possible  the  first  thorough  analysis  of  interference  effects  in  synthetic  aperture  im¬ 
ages  of  closely  spaced  scatterers.  Effects  of  errors  in  motion  compensation  are  also 
discussed  in  considerably  greater  depth  than  in  earlier  work.  Clomputer  simulation 
has  provided  us  with  a  useful  tool  with  which  to  verify  our  analytical  results;  we 
describe  our  simulation  method  and  present  a  comprehensive  selection  of  simulated 
images  illustrating  and  confirming  our  theoretical  analysis.  Finally  some  practical 
implications  of  our  results  for  the  autofocusing  of  ISAR  images  are  discussed. 

This  work  was  carried  out  under  Defence  Research  Package  AS02BM23 
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1  Introduction 


In  a  recent  memorandum  [l]  we  described  a  unified  theory  of  SAR,ISAR 
and  hybrid  SAR/ISAR  imaging  of  targets  which  took  account  of  target  rotational 
motions  and  the  relative  centre  of  mass  motions  of  the  target  and  the  radar  plat¬ 
form.  Some  preliminary  simulations  of  the  ISAR  imaging  of  a  point  target  were  also 
presented  and  were  found  to  be  in  excellent  agreement  with  the  predictions  of  our 
formal  analysis. 

In  this  paper  we  will  discuss  the  implications  of  this  work  for  ISAR  imaging 
in  considerable  detail.  After  reviewing  the  basic  concepts  and  formalism  of  [1]  we 
apply  these  methods  to  analyse  the  principles  of  ISAR  more  clearly  and  in  much 
greater  depth  than  has  been  done  in  pevious  work.  Particular  attention  is  paid  to  the 
properties  of  the  complex  (I,Q),  rather  than  the  real,  intensity,  image  so  that  we  are 
able  to  describe  coherent  interference  effects  in  images  of  closely  spaced  scatterers 
for  the  first  time.  Image  degradation  due  to  errors  in  motion  compensation  is  also 
discussed  exhaustively.  In  both  these  cases  the  powerful  techniques  introduced  in 
[1]  once  again  provide  us  with  a  concise,  physically  transparent  description  of  effects 
which  would  be  prohibitively  difficult  to  analyse  within  the  context  of  standard  ISAR 
theory. 

Computer  simulation  provides  us  with  a  useful  tool  with  which  to  confirm 
our  theoretical  findings.  Consequently  we  describe  our  methodology  in  more  detail 
than  was  done  in  [1]  and  present  an  extensive  set  of  simulated  images  demonstrating, 
and  confirming  in  detail,  our  theoretical  results.  Finally  we  discuss  the  practical 
implications  of  our  findings  for  the  autofocusing  of  ISAR  images. 


2  The  k-space  description  of  synthetic  aperture 
imaging 


In  our  earlier  work  [1]  we  presented  a  general  theoretical  framework  within 
which  we  can  describe  synthetic  aperture  imaging  concisely  and  in  detail.  The 
application  of  these  methods  to  the  special  case  of  inverse  synthetic  aperture  radar 
(ISAR)  is  discussed  in  this  section,  permitting  us  to  simplify  and  clarify  several 
features  of  our  original  approach.  Our  initial  discussion  also  serves  to  make  the 
present  work  essentially  self  contained,  by  providing  background  material  for  our 
subsequent  discussions  of  interference  effects  and  the  defocusing  and  distortion  of 
images  arising  from  errors  in  motion  compensation  and  allowing  us  to  discuss  our 
detailed  simulation  studies  in  the  overall  context  of  the  k-space  formalism. 
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2.1  Evaluation  of  the  Point  Spread  Function 


An  appealing  feature  of  the  k-space  formalism  we  apply  to  the  imaging 
process  is  the  way  in  which  its  component  parts  correspond  directly  to  operations 
carried  out  in  the  process  itself.  Thus  in  the  case  of  ISAR  a  signal  propagates  at 
time  t  with  wavenumber  |  (the  factor  of  |  is  included  for  subsequent  convenience) 
from  a  source  at  a  distance  r(<)  along  a  line  of  sight  vector  f  to  a  point  target  at 
r'(t)  ,  whose  motion  is  due  to  its  rotation  about  the  origin  of  our  coordinate  system. 
In  a  monostatic  configuration  the  signal  propagates  back  to  r(t)f  and  is  received  as 

Si{k,t)  =  exp(zfc|r(<)f-  r'(t)l)  (1) 

The  rotational  motion  of  r'(<)  can  be  represented  by 

r'(f)  =  M.(0.r'  (2) 


where  r'  is  the  position  of  the  imaged  point  at  time  zero.  M,(<)  is  the  linear  operator 
generating  the  target  rotation  in  time  t  [1,2].  As  this  is  orthogoneil  we  have 


\r{t)T -T'{t)\  =  |r(t)-r' 


(3) 


where 

Thus  we  may  write 


T{t)  =  r(t)M:\t).T 

Si{k,t)  =  exp(zfc|r(t)  —  r'|). 


(4) 

(5) 


(The  return  from  an  extended  target  may  be  derived  by  integrating  r'  over  the  spatial 
extent  of  the  target,  with  a  suitable  scattering  function  weighting  [1].)  Knowing  the 
range  r{t)  we  can  subject  this  signal  to  motion  compensation  to  give 


S2{k,t)  =  exp(-ifcr(f)  +  iA:|r(t)  -  r'j) 


(6) 


If,  as  in  [1],  we  consider  imaging  in  the  far  field,  we  may  approximate  this  by 

52(A’,f )  =  exp(— ik(<).r')  (7) 

where  k(t)  =  A:M7*(t).r,  so  that  returns  at  successive  times  t  provide  a  set  of  spatial 
Fourier  components  of  the  point  target’s  position  at  time  zero,  with  wave  vectors 
selected  by  its  own  rotational  motion  and  by  the  frequency  modulation  applied  to 
the  transmitted  waveform.  Thus,  if  the  rotation  of  the  target  is  specified  by  an  axis 
n,  and  an  angular  velocity  u:  [1,2]  we  have 


k(<)  =  (fco  +  ^k)  ((1  -  cos(u;<))nn.f  +  cos(wt)r  +  sin(a;t)(n  A  f)) . 


Here 


ko  = 


47r/o 

c 


(8) 
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where  fo  is  the  centre  transmitted  frequency  and  c  is  the  velocity  of  light,  6k  is  the 
change  in  transmitted  wave  number  due  to  frequency  modulation  and  n  =  — n,  (c.f. 
(4)).  On  expanding  (8)  to  second  order  in  6k  and  we  find  that 


k(t)  =  koT  +  koit’tn  A  f  -t-  6ku>tn  A  r  -f  A  (n  A  f )  (9) 

Consequently  k  space  is  interrogated  on  the  surface  of  a  cone  with  axis  n  and  con¬ 
taining  ko  [1].  In  [1]  we  showed  that  conventional  ISAR  processing  is  equivalent  to 
the  Fourier  inversion  of  these  components,  assuming  them  to  be  sampled  uniformly 
on  a  Cartesian  grid  whose  1,2  axes  (©1,62)  are  parallel  to  f, n  Ar  respectively.  (This 
coordinate  sytem,  supplemented  by  a  third  orthogonal  axis  63  parallel  to  f  A  (n  A  f ) 
will  be  adopted  throughout  our  subsequent  discussion.)  We  now  make  the  following 
identifications: 

ko  =  koT  (10) 

1  r  1  -  1  ^0  / 1  *1  \ 


k^i  =  6kr  =  kni  ■ 


k«2  =  A  ko  =  kfj2 


n  A  ko 


2  ^  ^  ^  J  1  •  Ol  ' - ' 

ko\  Sin  ll\ 

where  fl  is  the  angle  between  n,ko.  One  dimensional  Fourier  inversion  with  respect 
to  kffi  is  equivalent  to  the  acquisition  of  range  resolution  by  pulse  compression  while 
subsequent  inversion  with  respect  to  kH2  is  equivalent  to  the  Doppler  component  of 
ISAR  processing.  6k  and  t  may  be  eliminated  from  the  remaining  terms  in  (9)  to 
give 

k{t}  =  ko  +  kH  +  6k{kH)  (13) 

where  ^  ^ 

^k(kg)  =  ^ A  (n  A  r)  (14) 

ko  2  ko  sin  U 

^  kHikH2 . ^  ^  1  (cos  n  sin  Dca  -  sin^  Dei)  (15) 

ko  2  ko  sin  D  '  ' 

(This  straightforward  analysis  yields  the  result  expressed  in  [l]  as 


6ka  —  ^^'a^fkffQkfJ^^ 

and  allows  us  to  identify  elements  of  the  curvature  tensor  Ca0i  immediately.  Here 
and  elsewhere  we  adopt  the  convention  of  summation  over  repeated  indices.) 

By  bringing  these  results  together  we  may  write  the  point  spread  function 
of  the  ISAR  system  as 

H{rH,T')  =  (16) 
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H{'kji)  is  the  weighting  in  k-space  resulting  from  the  system’s  beam  and  bandwidth 
and  is  modelled  conveniently  by  the  Gaussian  function 


iy(kH)  = 


exp  (-iicjf.D.kff) 
(27r)(detD-')» 


where 


dll  0 
0  dia 


The  phase  term  exp(— tr'.6k)  can  be  expressed  similarly  as  exp  |k^.C.k/f)  where 


Cll  Ci2 
C21  C22 


with  C0y  =  r'^Ca0y  =  Cy0  or,  more  explicitly 


12. 

kfi 


Thus  we  may  express  the  point  spread  function  in  the  form 

f  (-5k*.(D  +  iCj.kfl)  (20) 

(27r)3(det  ^  ' 

g«ko.ArH  _ 27r _ r  7  1  \ 

~  (27r)2(detD-i)^)  [(q  -  /3)2  +  V  2“ 

X  exp  ^-^Ar^.A.Ar/rj  exp  ^--Afn-P-Ar/f  j  (21) 

where  a  =  ^11^22,  0  =  C11C22  —  cjj,  7  =  diiC22  +  d22Cn  and  At//  —  Th  —  r'jj  , 


(a  —  /?)AdjD  +  7AdjC 
(a  -  0y  +  7* 


Ia~0)d2i+-rc2j  -‘fC2i 

{a-0)^+ll^  (a-/?P+7’ 

-■ycia  {a-0)du+'t^ii 

{a-0)^+-ri  {a-0)^+y^ 
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and 


(a  —  yS)AdjC  +  7AdjD 
(q  -  I3y  +  72 


ia-0)ci2—f<ii3 

-(a-/?)cn 


~{a-0)c2i 

(a-/3)»+7» 

(a-0)cii  — yJn 
(<.-/9)»+7* 


The  amplitude  variation  within  the  point  spread  function  is  contmned  in  the 
Gaussian  function  exp  j  while  the  phase  variation  is  contained  in 

the  term  exp  |Ar^.P.Arfl^) .  In  [1]  we  have  shown  how  the  lower  order  moments 
of  the  point  spread  function  are  related  simply  to  the  elements  of  C  and  D  through 

((Aa:)=)  =  ^  (24) 


{AxAy)  = 


C12C22 


*****'^  2  2 

{{Ayf)  =  d22  +  ^  ^  (26) 

“22  “11 

The  confirmation  of  these  results  by  computer  simulation  is  presented  in  [l]. 
The  compact  representation  of  the  phase  variation  within  the  point  spread  function 
enables  us  to  describe  interference  effects  within  synthetic  aperture  images  of  closely 
spaced  scatterers  by  the  appropriate  coherent  summation  of  point  spread  functions 
centred  on  the  positions  of  the  scatterers.  The  amplitude  and  phase  of  the  resulting 
function  will  then  characterise  the  detailed  structure  of  the  corresponding  ISAR 
image.  This  procedure  will  be  illustrated  by  computer  simulation  in  a  subsequent 
section  of  this  Memorandum. 


2.2  Motion  Compensation  Errors 


The  C  matrix  formalism  we  have  just  described  is  sufficiently  flexible  to 
accommodate  defocusing  effects  due  to  errors  in  motion  compensation  i.e.  in  the 
estimation  of  the  position  of  the  centre  of  rotation  of  the  target.  A  preliminary  dis¬ 
cussion  of  this  general  problem  was  given  in  [1];  by  restricting  our  present  discussion 
to  inverse  synthetic  aperture  radar  we  are  able  to  highlight  several  points  which  were 
neglected  in  our  earlier  work.  Once  again  our  formalism  has  become  more  stream¬ 
lined  so  that  much  of  the  algebraic  drudgery  that  was  evident  in  [l]  can  be  avoided. 


Let  us  assume  that  we  have  made  an  error  6r(/)  in  estimating  the  position  of  the 
centre  of  rotation  of  the  target.  The  returned  signal  Si{k,t)  is  now  given  by 

Si{k,t)  =  exp{ik\r{t)T  -  M,(t).r'  -  6r(t)|).  (27) 

which,  when  subjected  to  motion  compensation  based  on  the  range  estimate  r{i) 
yields  (in  the  far  field  limit) 

S2ik,t)  =  exp(— tk(t).r' —  ik(t).M7\i)'^K0)  (2®) 

Thus  we  see  how  the  phase  of  the  processed  signal  has  been  altered  as  a  result  of 
incorrect  motion  compensation.  This  phase  change  results  in  a  shift  in  the  posi¬ 
tion  of  the  centre  of  the  point  spread  function,  which  is  also  subject  to  a  further 
defocusing.  Quantitative  measures  of  these  effects  can  be  obtained  by  expanding 
this  phase  change  to  0{kjj)]  the  linear  term  will  specify  the  image  shift  while  the 
quadratic  terms  characterise  the  extra  defocusing  and  may  be  represented  concisely 
by  a  suitable  augmentation  of  the  elements  of  the  C  matrix  (19).  As  before  we  write 

k(<)  =  (1  -I-  ^)  ^ko  +  u;t(n  A  ko)  -I-  A  (n  A  ko)j  .  (29) 

If  6r(t)  is  expressed  in  terms  of  a  gross  displacement  error  ^r,  a  velocity  v  and  an 
acceleration  a  through 

6r(<)  =  +  vf -f -at*  (30) 

we  find  by  direct  expansion  that 

M7^(t).^r(t)  =  -f  (v  -(-  n  A  ^r)t  +  -(a  -f  2u;n  A  v  -|-  u;*n  A  (n  A  6r))t*  (31) 

where  v  and  a  have  been  supplemented  by  terms  corresponding  to  the  centrifugal 
and  Coriolis  corrections  applied  in  the  discussion  of  dynamics  referred  to  a  rotating 
coordinate  frame[3].  By  expanding 

k{t)M:\t)Mt)  =  (1  +  ^)  ^ko  +  u;t(n  A  ko)  +  A  (n  A  ko)j 

.  (v  -I-  n  A  6r)t  -f  ^(a  -f  2u;n  A  v  +  u7*n  A  (n  A  6r))t*  j  ,  (32) 

exploiting  the  vector  identities 

(n  A  (n  A  ko).^r)  =  (n  A  (n  A  6"r).ko  =  -(n  A  ko).(n  A  ^r) 
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and  eliminating  t  (c.f.  (12))  we  find  that 


k(t).M,  ^(t).^r(t)  =  ko.^r  -t-  kniSrHi 

ko.vfeg2  (ko.v)fejyifejy2  1  kp.afejf; 

u;|n  A  kol  fcott>|nAko|  2u?2|nAkoP’ 

noting  the  exact  cancellation  of  terms  arising  from  the  centrifugal  and  Coriolis  cor¬ 
rections.  The  point  spread  function  incorporating  the  effects  of  errors  in  motion 
compensation  is  given  by 


(27r)^  J 


(34) 


On  comparing  this  result  with  (16)  we  see  that  the  presence  of  the  zeroth  order 
term  in  (33)  merely  multiplies  the  point  spread  function  by  a  constant  phase  factor 
exp(— iko-^r).  The  linear  term  in  (33)  shifts  the  centre  of  the  point  spread  function 
from 

^  H 


to 


Th  =  r'^  -f  btHiii  + 


re2 


u;|  sin  fl| 

while  the  quadratic  terms  modify  the  defocusing  to  that  represented  by  the  C  matrix 


®  fco  +  u;|s”inn|) 


(35) 


Thus  we  see  that  the  position  of  the  image  is  shifted  by  the  along  range  part  of  the 
gross  error  in  the  estimation  of  the  target /platform  separation,  and  across  the  line  of 
sight  by  the  uncompensated  along  range  velocity.  This  observation  is  consistent  with 
the  ISAR  technique  acquiring  across  range  resolution  only  through  the  rotational 
motion  of  the  target  and  so  being  incapable  of  detecting  the  constant  6r2.  We  also 
note  that  the  modified  C  matrix  (35)  is  independent  of  ^r.  This  is  a  physically 
reasonable  result  as  we  can  see  by  considering  the  special  case  of  a  point  target 
at  the  origin  of  the  coordinate  system.  Merely  misplacing  such  a  target  (which  is 
not  affected  by  rotation)  cannot  cause  defocusing;  our  inclusion  of  centrifugal  and 
Coriolis  type  terms  in  the  analysis  has  eliminated  the  dependence  of  C  on  8t  and 
so  has  ensured  that  our  formal  and  intuitive  understandings  of  the  imaging  process 
are  in  accord. 
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From  (24-26)  we  see  that  the  lower  order  moments  of  the  point  spread 
function,  corrupted  by  incorrect  motion  compensation,  are  given  by 


1 


((Ai)  )  =  dll  -4-  ,  ,2  I  7“  ~T~- — o 

a22«o  \«o  u>|smf2 


Vl 


{Ax  Ay)  =  -- 


Vl 


^22^0  \A:o  u;|sinfi| 


rj  —  Tj  cot  U  — 


Cl 


u;*  sin^ 


fi) 


(36) 

(37) 

■) 

(38) 

These  results  furnish  us  with  a  quantitative  measure  of  the  defocusing  associated 
with  incorrect  motion  compensation;  in  particular  we  note  the  results  which  de¬ 
scribe  defocussing  due  to  Fourier  component  sampling  on  a  conical  surface  are  not 
supplemented  additively  but  in  a  more  subtle  fashion,e.g. 


((Aj,)»)  =  1  (j-  (r;  -  r^cot  ft  -  ^ 


Vl 


u}\  sin 


((Ax)2)^dii  + 


+ 


^22^0  ^22^0^^  sin^  n 


The  extension  of  standard  ISAR  theory  to  include  these  effects  would  be  difficult, 
both  conceptually  and  in  terms  of  its  algebraic  complexity. 


3  Simulation  Results 


In  this  section  we  describe  computer  simulations  which  verify  the  theoretical 
results  discussed  in  the  previous  section.  Simulated  images  of  a  ship  model  illustrate 
how  different  target  rotations  generate  images  projected  onto  different  planes.  The 
variation  of  the  phase  of  the  point  spread  function  is  ansilysed  in  detail,  allowing  us 
to  describe  the  interference  effects  in  images  of  closely  spaced  scatterers  observed 
in  simulated  ship  images.  Finally  we  demonstrate  the  effects  of  errors  in  motion 
compensation  on  the  amplitude  and  phase  of  the  point  spread  function  and  the 
more  pronounced  interference  effects  observed  in  defocused  images.  Throughout  the 
agreement  we  obtain  between  theory  and  simulation  is  excellent  and  provides  a  very 
satisfactory  confirmation  of  our  formal  analysis. 


3.1  Simulation  Methods 

The  simulation  procedure  we  adopt  consists  of  two  stages:  the  generation  of  a  succes¬ 
sion  of  target /platform  configurations  and  the  evaluation  of  the  monostatic  returns 
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in  each  of  these  configurations  for  a  range  of  transmitted  frequencies,  and  their  sub¬ 
sequent  processing  to  yield  a  simulated  ISAR  image.  As  the  successive  steps  in  the 
imaging  process  are  mirrored  very  closely  in  our  formal  development  we  will,  where 
appropriate,  employ  a  notation  based  on  that  of  the  previous  section.  Through¬ 
out  we  assume  that  the  radar  platform  is  stationary  (i.e  that  r{t)  =  r);  referred  to 
our  standard  coordinate  system  it  has  a  position  vector  r  =  (r,0,0).  The  target 
consists  in  general  of  a  set  of  point  scatterers,  labeled  by  index  m.  To  each  point 
scatterer  we  assign  an  initial  position  its  subsequent  positions  are  generated 

by  the  straightforward  numerical  implementation  of  (2), (30)  i.e.  at  time  ti  =  lAt 
the  position  of  the  scatterer  is  given  by 

(6r  +  yti  -f  •+■  M,(</).r'„ 

M,(<j)  generates  the  target  rotation  through 

=  (1  -  cos(wt/))n,(n,.rm)  -I-  cos(wt/)r',„  -f  sin(u;</)n,  A  (39) 

At  each  time  ti  a  set  of  motion  compensated  returns  is  calculated  as 

=  exp(-ikjr)  €xp{ikj\T  -  St{U)  -  Mg(t/).r'„|)  (40) 


the  wavenumber  kj  being  given  by 


—  ko  -\- 


47rjA/ 

c 


where  A/  is  the  frequency  modulation  increment.  By  stepping  through  equal  num¬ 
bers  of  positive  and  negative  integer  values  of  j,  I  and  summing  over  all  the  point 
scatterers  in  the  target  a  256  x  256  array  of  backscattered  returns  is  built  up.  This 
is  then  weighted  with  a  Gaussian  function  centered  on  ji  =  /  =  0,  corresponding  to 
the  weighting  function  Hiku).  To  simulate  the  imaging  process  we  employ  a  two  di¬ 
mensional  Fast  Fourier  transform  to  effect  the  inversion  of  this  weighted  array.  The 
equivalence  of  this  procedure  to  standard  ISAR  processing  [4]  is  discussed  in  [l];  this 
relationship  is  represented  schematically  in  Figure  1.  The  frequency  modulation,  ex¬ 
perimental  geometry  and  target  rotation  rate  define  the  length  scales  characteristic 
of  the  image,  the  1  coordinate  being  discretised  in  units  of 


xstep  = 


c 

512A/ 


and  the  2  coordinate  in  units  of 

c 

^  512Af/ou;,|sinf2l 


(In  comparing  simulation  and  theoretical  results  care  must  be  taken  to  ensure  that 
distances  within  the  image  plane  are  always  scaled  appropriately  in  these  units  and 
that  the  convention  n  =  — n,  is  observed  throughout.) 
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3.2  Ship  Images 


Our  first  set  of  simulation  results  illustrates  the  way  in  which  different  target 
rotational  motions  or,  more  specifically,  the  axes  about  which  these  rotations  occur, 
generate  ISAR  images  projected  onto  different  planes.  The  axis  of  image  projection 
can  be  thought  of  as  that  axis  along  which  we  have  no  information  about  position. 
We  recall  that  the  range  profile  provides  us  with  information  along  r,  while  the 
Doppler  profile  provides  us  with  information  along  (n,  A  r).  Consequently  the  axis 
of  projection  must  be  normal  to  both  these  vectors  and  so  lies  along  r  A  (n,  A  r)  i.e. 
the  image  is  always  formed  in  the  (1,2)  plane  of  our  standard  coordinate  system.  We 
have  simulated  images  of  a  250  scatterer  model  of  an  F88  frigate  undergoing  different 
rotational  motions,  which  we  present  along  with  line  drawings  of  the  ship  model, 
projected  onto  the  (1, 2)  planes  defined  by  these  motions.  These  are  shown  in  Figures 
3-6,  where  yaw,  pitch  and  composite  rotations  can  be  seen  to  produce  plan,  side,  45® 
elevated  side,  and  three  quarters  views  of  the  ship  .  Table  1  summarises  the  radar 
parameters  and  ship  motions  used  to  produce  each  image,  allowing  a  comparison 
between  their  theoretical  and  actual  planes  of  projection.  Visual  inspection  reveals 
that  these  planes  are  indeed  coincident. 


3.3  Phase  Effects 

A  more  exacting,  quantitative  test  of  our  theory  is  provided  by  the  com¬ 
parison  of  our  point  spread  function  results  (21-23)  with  the  simulated  images  of  a 
single  point  scatterer.  In  our  earlier  work  [1]  we  demonstrated  the  satisfactory  agree¬ 
ment  obtained  between  theory  and  simulation  for  the  amplitude  variation  within  the 
image  of  a  single  point.  We  now  demonstrate  that  the  phase  variation  within  such 
images  is  also  well  described  by  these  results.  In  Figures  8,9  we  present  a  compar¬ 
ison  of  the  variation  in  phase  in  simulated  images  obtained  as  described  in  Table  2 
and  our  corresponding  theoretical  predictions.  (Throughout  this  report  we  present 
our  results  in  a  colour-coded  contour  plot  format.  Amplitude  will  be  displayed  in 
arbitrary  units  (which  will,  of  course,  be  the  same  in  comparing  any  particular  set 
of  simulated  and  theoretical  results)  while  phase  will  be  presented  in  radians  and  be 
taken  to  lie  in  the  range  [— 7r,-|-7r].  We  also  draw  the  reader’s  attention  to  the  fact 
that,  in  Table  2,  scatterer  position  vectors  are  presented  in  Cartesian  component 
form;  all  other  vectors  are  presented  the  (r,  6,  <f))  spherical  polar  form  illustrated  in 
Figure  2.)  The  agreement  can  be  seen  to  be  very  satisfactory.  It  should  be  stressed 
that,  in  obtaining  these  results,  care  was  taken  to  ensure  that  the  simulations  were 
carried  out  in  the  far  field  limit,  i.e.  that  (if  is  the  component  of  the  target’s  posi¬ 
tion  vector  perpendicular  to  the  line  of  sight  vector  (1,0,0))  the  Fraunhofer  criterion 


was  satisfied. 


In  the  course  of  our  ship  image  simulations  it  was  noticed  that  when  a 
target  contained  several  point  scatterers  placed  in  close  proximity  the  image  obtained 
frequently  showed  features  most  readily  interpreted  as  interference  effects.  We  now 
present  results  analysing  this  phenomenon  in  more  detail.  If  we  have  more  than 
one  scatterer  in  the  target  we  will  obtain  an  image  given  by  the  coherent  sum  of 
the  complex  point  spread  functions  (21)  centred  on  the  positions  of  these  point 
scatterers,  which  can  in  general  display  interference  effects  manifested  as  intensity 
maxima  and  minima.  Thus  the  form  of  the  complex  image  of  a  set  of  scatterers, 
labelled  by  index  m,  is  given  by 

/(■•»)' =  («) 

m 

where  (c.f.  (21))  each  of  the  depends  on  the  position  vector  r'^  of  the 

associated  point  scatterer.  The  complicated  dependence  of  the  phase  and  amplitude 
of  the  resulting  complex  image  on  the  target  and  experimental  geometry  can  lead  to 
quite  subtle  features  appearing  in  the  imagery,  differing  from  those  expected  from  a 
superposition  of  the  intensities  of  the  signatures  of  the  individual  point  scatterers. 
This  effect  is  illustrated  in  Figure  10.  Figure  10a  shows  the  intensity  image  of  a 
point  scatterer  at  (0,100,0)  relative  to  the  centre  of  mass  of  the  target,  rotating 
about  an  axis  along  the  3-axis  while  Figure  10b  shows  the  simulated  intensity  image 
of  two  points  at  (0,100,0)  and  (0,101,0).  Figure  10c  show  the  results  obtained  by 
adding  the  intensities  of  images  of  points  at  (0,100,0)  and  (0,101,0)  incoherently. 
This  does  not  agree  with  the  simulation.  However  the  theoretical  result  (Figure 
lOd)  based  on  Equation  (42)  is  found  to  be  in  close  agreement  with  the  simulation. 
It  was  also  noted  that  the  simulated  interference  pattern  becomes  asymmetric  when 
the  radar  is  taken  out  of  the  far-held  limit  i.e.  when  the  inequality  (41)  was  not 
satisfied.  These  near-held  effects,  resulting  from  the  curvature  of  the  wave  front 
incident  on  the  target,  have  recently  been  analysed  in  considerable  detail  and  are 
discussed  elsewhere  [6]. 

Finally  Figures  11  and  12  show  the  excellent  agreement  between  simulation 
and  theory  for  ten  closely  spaced  scatterers.  Figures  11a  and  11b  illustrate  the 
agreement  between  the  amplitude  variation  whilst  Figures  12a  and  12b  show  the 
phase  variation.  Taken  as  a  whole  the  results  we  have  presented  so  far  provide 
adequate  conhmation  of  the  validity  of  our  description  of  the  point  spread  function 
of  an  ISAR  system  for  which  motion  compensation  is  perfect. 
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3.4  Effects  of  Incorrect  Motion  Compensation 


Our  final  set  of  simulation  results  illustrates  and  confirms  bur  analysis 
of  the  effects  of  incorrect  motion  compensation.  Figure  13  shows  the  amplitude 
distribution  of  the  image  of  a  non-rotating  point  target  with  an  uncompensated 
velocity  and  acceleration.  Once  again  colour-coded  contour  plots  of  the  simulated 
images  are  presented,  along  with  plots  of  the  appropriate  theoretical  results.  The 
corresponding  phase  plots  are  shown  in  Figure  14.  These  results  can  be  compared 
with  those  presented  in  Figures  15  and  16,  which  show  the  amplitude  and  phase 
plots  for  the  image  of  a  rotating  point  target  with  the  same  uncompensated  velocity 
and  acceleration.  The  excellent  agreement  between  theory  and  simulation  is  evident 
both  visually  and  in  the  numerical  results  given  in  Table  3  for  (Ax^),  (Ay^),  and 
{Ax Ay).  The  differences  between  Figures  13,14  and  15,16  illustrate  the  effect  of 
the  coupling  between  ‘conical  polar  formatting’  and  centre  of  mass  motion  effects 
manifest  in  the  results  (36-38). 

Figure  17a  shows  the  amplitude  distribution  plot  for  the  image  of  two  closely 
spaced  point  targets  assuming  perfect  motion  compensation.  Theory  (17b)  and  sim¬ 
ulation  are  presented  and  can  be  seen  to  be  in  agreement.  The  corresponding  plots 
for  the  image  of  equivalent  point  targets  rotating  in  a  coordinate  system  with  an 
uncompensated  velocity  and  acceleration  are  shown  in  Figure  17c/d.  The  effects  of 
centre  of  mass  motion  in  defocussing  the  image  and  introducing  discernible  interfer¬ 
ence  features  are  evident  in  the  simulated  results,  which  can  be  seen  to  be  described 
accurately  by  the  theory.  Table  2  shows  the  parameters  used  in  these  simulations. 


4  Autofocusing 


The  ultimate  objective  of  this  work  is  to  establish  a  capabilty  to  make  better  use  of 
radar  data  for  target  localisation  and  identification.  It  is  therefore  desirable  to  derive 
well  focused  images  in  situations  where  motion  parameters  are  only  known  approx¬ 
imately.  One  approach  to  this  is  to  estimate  the  unknown  motions  by  maximising 
the  focusing  of  the  image  with  respect  to  each  of  the  unknowns.  Thus  for  SAR 
the  across-track  acceleration  is  often  estimated  by  maximising  the  image  contrast  or 
minimsing  tjie  sub-aperture  spatial  offset.  For  more  complicated  problems,  such  as 
imaging  a  moving  and  rotating  ship  from  an  airborne  sideways  looking  radar,  the  un¬ 
known  parameter  space  has  too  many  dimensions  and  subtle  defocusing  interactions 
between  motions,  to  allow  a  practical  image  optimisation. 

The  analysis  presented  in  earlier  sections  of  this  report  and  [1]  shows  that  the 
defocusing  of  SAR,  ISAR  and  hybrid  SAR/ISAR  images  can  be  expressed  concisely 
in  a  matrix  form.  Manipulation  of  the  results  show  how  the  following  effects  cause 
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image  distortion  and  defocusing, 

(a)  Approximations  in  the  processing,  e.g.  ‘Range-Doppler  processing  , 

(b)  Unwanted  motions,  e.g.  rotational  accelerations, 

(c)  Errors  in  estimations  of  other  motion  parameters,  e.g.  target  rotation 
rate,  target  and  radar  velocities  and  accelerations,  and 

(d)  Scatterers  on  the  target  offset  from  the  imaging  plane. 

A  particular  strength  of  the  approach  is  the  grouping  together  of  terms 
which  cause  similar  effects  on  imagery,  and  are  therefore  ambiguous  to  interpretation. 
Simple  examples  of  this  are 

(a)  Rotation  rate  and  image  scaling  in  ISAR,  and 

(b)  The  first  order  effect  of  line-of-sight  motion  and  azimuthal  position  in 
SAR  imagery. 

The  formalism  therefore  allows  the  search  for  autofocusing  to  be  undertaken 
within  the  correct  unambiguous  parameter  space.  Also,  specific  functional  relation¬ 
ships,  between  the  degree  of  defocusing  and  the  position  on  the  target,  are  associated 
with  each  type  of  motion.  This  suggests  that  some  decoupling  of  searches  of  sets 
of  parameters  may  be  possible.  Finally  the  quantitative  linkage  between  defocusing 
and  actual  motion  should  provide  bounds  within  which  to  undertake  the  searches, 
and  measures  of  the  accuracy  of  the  focus  achieved. 

Specific  algorithms  have  not  yet  been  proposed,  nor  have  focus  optimisa¬ 
tion  measures  been  tested.  One  particular  concern  relates  to  the  use  of  contrast 
maximisation,  since  one  finding  of  this  work  has  been  the  large  phase  variations 
associated  with  defocused  point  spread  functions.  This  may  result  in  high  contrast 
fringes  between  large  defocused  scatterers  of  the  type  we  have  analysed  in  section 
3.4.  Future  work  will  address  these  features,  with  the  aim  of  developing  a  well 
founded  SAR/ISAR  autofocusing  method,  which  performs  better  than  the  current 
ad~hoc  approaches  and  provides  a  higher  yield  of  classifiable  target  images. 


5  Conclusions 


In  this  Memorandum  we  have  discussed  many  features  of  far  field  ISAR  images  in 
much  greater  detail  than  has  been  done  previously.  This  has  been  made  possible  by 
the  application  of  the  general  theory  of  synthetic  aperture  imaging  introduced  in  our 
earlier  work  [1].  Our  present  work  has  confirmed  the  utility  and  power  of  this  method 
and  at  the  same  time  has  simplified  and  clarified  many  of  the  calculational  procedures 
involved  in  its  application  to  specific  problems.  Extensive  computer  simulations  have 
confirmed  its  value  in  the  analysis  of  the  phase  characteristics  of  the  point  spread 
function  and  the  occurence  of  interference  effects  in  synthetic  aperture  images  of 
closely  spaced  scatterers.  An  essentially  complete  description  of  the  effects  of  errors 
in  motion  compensation  in  far-held  ISAR  imaging  has  also  been  achieved,  and  has 
again  been  verified  by  computer  simulation.  These  developments  have  furnished  a 
very  sound  foundation  for  subsequent  work  in  which  ISAR  images  will  be  enhanced 
by  autofocus  methods.  Consequently  we  have  presented  a  preliminary  assessment  of 
the  difficulties  likely  to  be  encountered  in  implementing  these  techniques  and  their 
evident  amelioration  by  the  results  of  our  work  to  date. 
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Table  Of  Figures 

Figure  1  A  diagrammatic  summary  of  the  simulated  imaging  process. 

Figure  2  The  rotational  coordinate  system  used  in  the  theory  and  the  simulations. 

Figure  3  A  simulated  image  of  a  point  scatterer  modelled  ship  in  comparison  with 
a  plan  view  of  the  ship.  In  this  case  the  ship  is  yawing  at  a  rate  of  j sec. 

Figure  4  A  simulated  image  of  a  point  scatterer  modelled  ship  in  comparison  with 
a  side  view  of  the  ship.  In  this  case  the  ship  is  pitching  at  a  rate  of  4  j sec. 

Figure  5  A  simulated  image  of  a  point  scatterer  modelled  ship  in  comparison  with 
a  45“  elevated  side  view  of  the  ship.  The  ship  is  observed  along  the  x-axis  and 
its  motion  combines  yaw,  pitch,  and  roll. 

Figure  6  A  simulated  image  of  a  point  scatterer  modelled  ship  in  comparison  with 
a  three-quarters  view  of  the  ship.  The  ship  is  observed  along  the  —  xj/-axis 
and  its  motion  combines  yaw,  pitch,  and  roll. 

Figure  7  Rotation  axis  of  simulation. 

Figure  8  Phase  plot  of  image  of  scatterer  at  (67,67,0). 

Figure  9  Simulated  and  theoretical  images  for  two  point  scatterers. 

(a),(b)  Point  at  (67,67,0),n=(4,0,45).  (c),(d)  Point  at  (134,67, 134),n=(4,0,60) 

Figure  10  The  effects  of  coherent  and  incoherent  addition  of  scatterers,  and  com¬ 
parison  with  simulation. 

Figure  11  Amplitude  plot  showing  the  interference  between  10  scatterers.  Details 
in  Table  2. 

Figure  12  Phase  plot  showing  the  interference  between  10  scatterer  s.  Details  in 
Table  2. 

Figure  13  Amplitude  plot  showing  the  effect  of  two  centre  of  mass  motions  on  a 
scatterer  at  (0,0,0). 

Figure  14  Phase  plot  showing  the  effect  of  two  centre  of  mass  motions  on  a  scatterer 
at  (0,0,0). 

Figure  15  Amplitude  plot  showing  the  effect  of  two  different  centre  of  mass  motions 
on  a  scatterer  at  (67,67,0)  -  same  motions  as  Figure  13. 

Figure  16  Phase  plot  showing  the  effect  of  two  different  centre  of  mass  motions  on 
a  scatterer  at  (67,67,0)  -  same  motions  as  Figure  14. 

Figure  17  Interference  between  two  scatterers  (16,15,0)  and  (16,17,0),  with  and 
without  a  centre  of  mass  motion  of  the  target. 
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Table  1  -  Parameters  for  ship  images 


Figure  No 

Line  of  Sight 
k 

Rotation  axis 

n 

Theoretical 

Axis  — k  A  (n  A  k) 

Actual  Axis 
of  Projection 

3 

(1,0,0) 

(0,0,1) 

(0,0,-i) 

(0,0,-l) 

4 

(1,0,0) 

(0,1,0) 

(0,-l,0) 

(0,-l,0) 

5 

(1,0,0) 

(1,1,1) 

(0,-1, -1) 

(0,-Vi) 

6 

(1-1,0) 

(1,1,1) 

(-1,-1,-1) 

(-I,-l,-l) 

Rate  of  rotation  4  degs/sec. 

Centre  frequency  fo  =  lOGHz. 

Step  frequency  A  /  =  1  MHz. 

Time  Step  At=  1.8ms. 

S.d.  of  Gaussian  Weighting  25 

Motion  compensation  No  velocity  or  acceleration  correction. 
Radar  positioned  at  far  field  limit. 
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Table  2  -  Parameters  for  Figures 


Figure  No 

Rotation  axis 

n 

Displacement 

Ar 

Velocity 

V 

Acceleration 

a 

Position 

r' 

8 

(1,0,0) 

(0,0,0) 

(0,0,0) 

(0,0,0) 

(67,67,0) 

9a 

(1,45,0) 

(0,0,0) 

(0,0,0) 

(0,0,0) 

(67,67,67) 

9b 

(1,60,0) 

(0,0,0) 

(0,0,0) 

(0,0,0) 

(134,67,134) 

11 

(1,0,0) 

(0,0,0) 

(0,0,0) 

(0,0,0) 

(49.49.49)  (50,52,48) 

(53.50.50)  (51,51,0) 
(50,48,0)  (47,47,52) 
(52,48,30)  (50,50,30) 

(47.53.50)  (50,50,50) 

B 

B 

(1,34,50) 

(1,34,50) 

(1,57,63) 

(1,57,63) 

B 

(0,0,0) 

(67,67,0) 

m 

B 

B 

(4,214,63) 

(4,214,63) 

(2,70,60) 

(2,70,60) 

(0,0,0) 

(67,67,0) 

17a 

17c 

(1,0,0) 

(1,0,0) 

(0,0,0) 

(0,0,0) 

(0,0,0) 

(3,50,40) 

(16,15,0)  (16,17,0) 
(16,15,0)  (16,17,0) 

Rate  of  rotation  4  degs/sec. 

Coordinate  Systems  Position  vector  in  {x,y,z),  others  in  (r,9,<j>) 
Centre  frequency  /„  =  lOGHz. 

Step  frequency  A  /  =  3  MHz. 

Time  Step  A  t=  Sms. 

S.d.  of  Gaussian  Weighting  25 

Motion  compensation  No  velocity  or  acceleration  correction. 
Radar  positioned  at  far  field  limit. 


18 


Table  3  -  Moments  comparing  theory  and  simulation 


■ 

(0,0,0) 

FiKurelSa  FigurelSc 

(67,67,0) 

Figurel5a  Figuicl5c 

Sim. 

Theory 

1  Sim. 

1  Sim. 

1  Sim. 

Theory 

fl 

■ 

B 

■ 

fl 

■ 

B 

B 

0.393 

(y*) 

1.509 

1.511 

5.824 

5.835 

■ 

B 

4.426 

4.446 

(*y) 

■ 

fl 

■ 

B 

B 

B 

1.058 

1.062 

Simulation  Results  -  derived  numerically  from  image. 
Theoretical  Results  -  derived  from  equations  [44-46]. 
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Figure  2  The  rotational  coordinate  system  used  in  the  theory  and  the  simulations. 


Figure  3  A  simulated  image  of  a  point  scatterer  modelled  ship  in  comparison  with 
a  plan  view  of  the  ship.  In  this  case  the  ship  is  yawing  at  a  rate  of  4  /sec. 


Figure  4  A  simulated  image  of  a  point  scatterer  modelled  ship  in  comparison  with 
a  side  view  of  the  ship.  In  this  case  the  ship  is  pitching  at  a  rate  of  4°  j sec. 


Figure  5  A  simulated  image  of  a  point  scatterer  modelled  ship  in  comparison  with 
a  45"  elevated  side  view  of  the  ship.  The  ship  is  observed  along  the  s-axis  and 
its  motion  combines  yaw,  pitch,  and  roll. 


Figure  6  A  simulated  image  of  a  point  scatterer  modelled  ship  in  comparison  with 
a  three-quarters  view  of  the  ship.  The  ship  is  observed  along  the  — zy-axis 
and  its  motion  combines  yaw,  pitch,  and  roll. 


Figure  9  Simulated  and  theoretical  images  for  two  point  scatterers 

(a),(b)  Point  at  (67,67,0)m-(4,0.45).  (c).(d)  Point  at  (134,67,134),n  =  (4,0,G0) 


Figure  10  The  efietis  of  coherent  iiiul  incoherent  addition  ol  scatterers,  and  com¬ 
parison  with  simulation. 


Figure  13  Amplitude  plot  showing  the  eflect  ol  two  centre  of  mass  motions  on  a 
scatterer  at  (0,0,0). 


a]  Simulation 


c)  SimulBUon 


Figure  14  Phase  plot  showing  the  eilect  of  two  centre  of  mass  motions  on  a  scat 
tercr  at  (0.0.0). 


Figure  15  Amplitude  plot  showing  the  effect  of  two  different  centre  of  mass  motions 
on  a  scatterer  at  (67,67,0)  -  same  motions  as  Figure  13. 


a)  Simulation 


b)  Theory 


c)  Simulation 


Figure  16  Phase  plot  showing  the  eHect  ol  two  different  centre  of  mass  motions  on 
a  scatterer  at  (67.67,0)  -  same  motions  as  Figure  11. 


a)  Simulation  -  iKthout  cm  motion. 


b)  Theory  -  -without  cm  motion. 


lire  17  Interference  between  two  scatterers  (JtS. 1.5.0)  and  (16.J.7.0),  with  and 
without  a  centre  of  mass  motion  of  the  target. 
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